FINITE GELFAND PAIR APPROACHES FOR EHRENFEST DIFFUSION 

MODEL 
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Abstract. A classical diffusion model of Ehrenfest which consists of 2- urns and n-balls is 
realized by a finite Gelfand pair (H n , S n ), where H n is the hyperoctahedral group and S n is the 
symmetric group. This fact can be generalized to multi-urn version by using Gelfand pairs of 
complex reflection groups . 
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1. INTRODUCTION 

There are two urns, the left one containing n-balls and the right one having no ball. We 
shuffle the balls according to a rule as follows; At each step one ball is chosen randomly and 
the ball is moved to the other urn chosen randomly. This process is called Ehrenfest diffusion 
model 

We consider a generalization of the above process by increasing the number of urns. Let 
r > 2. There are r distinct urns Uq, U± 7 • • • , U r —\ : initially Uq contains distinct balls B\ ; • • • , B n . 
Let us consider a similar diffusion process as above: At each step, we choose a ball randomly 
and transfer the ball into one of the other urns. Here we consider some kind of destinations of 
picked balls. The first one is any other urns without where it was (Section 3). Another two 
cases are considered in the Appendix of this paper; The second case is that a ball picked from 
Ui is moved to U i+ i ( mod r \ (Appendix 14. The last case is that a ball picked from U{ is moved 
to U i+1 ( m odr) or Ui-.x ( mo d r) (Appendix SHJ) . 

If Uj contains the ball Bi, then we define a function on the balls by bi = b(Bi) = j. Then we 
can identify each configuration of our model with an element of a set 

B(r, n) = {(bi \ 1 < i < n) | < h < r - I}. 

In the next section, we see that the set B(r,n) is an realization of a finite homogenous space 
G(r,l,n)/S n , where G(r, l,n) is an complex reflection group and S n is a symmetric group. 
Indeed this pair of groups (G(r, l,n),S n ) is a Gelfand pair. Our main purpose is to analyze a 
stochastic space B(r,n) by using this Gelfand pair. Further the book [2] is good introduction 
for an application of finite Gelfand pairs to probability theory. 
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2. B(r,n) and (G(r,l,n),S n ) 

In this section we introduce how to identify B{r, n) with a certain finite and discrete homoge- 
nous space. Let S n be the symmetric group and G(r, l,n) = Z/rZ I S n the complex reflection 
group. We denote an element of G(r, l,n) by (x\,--- ,x n ;a), where Xi € Z/rZ and a € S n . 
Under this notation, we remark that S n is a subgroup {(0, • • • , 0; a) \ a G SV,,} of G(r, 1, re). We 
define an action of G(r, 1, n) on -B(r, n) by 

| 1 < i < re) = [xi + &o— (mod r) | 1 < z < re), 

where x = (a?i, • • • , x n ; a) € G(r, 1, n) and (6j | 1 < i < n) £ B(r, n). 

Proposition 2.1. This action is transitive on B(r,n). 

Proof. We set x = (r — bi, ■ ■ ■ , r — b n ;l). Then we have x(bi \ 1 < i < n) = (0, • • • , 0). 
Furthermore this action is invertible. Therefore any two elements are transferred to each other 
by the action of G{r, 1, n). □ 

Proposition 2.2. Put Iq = (0, • • • ,0) € B(r,n). Then the stabilizer of Iq is S n . 

Proof. It is clear that the definition of the action. □ 

From these propositions, we can identify B(r, n) with a finite homogenous space G{r, 1, n)/S n . 
Let N = {0, 1, 2, • • • }. For a composition m = (mo, • • • , m r _i) of n, let |m| = mo + • • • + m r _i 
be the size of m. Put N(r, n) = {m € N r | |m| = re}. For m € N(r, re), we define a partition by 
A(m) = (0 m °l mi • • • (r — l)" 1 "— 1 ), i.e, m« is regarded as the multiplicity of i. Let m\(xi, • • • , x n ) 
be a monomial symmetric polynomial indexed by a partition A. Let £ be a primitive rth root of 
unity. We set 

m A (£) = m A (i 1 _^i,e 1 _-^,--- , f-y- - L £3 

for (4,..- G iV(r,re). 

Then the following theorem holds. 

Theorem 2.3. [3] 

(1) A pair (G(r,l,n), S n ) is a Gelfand pair. 

(2) The permutation representation is decomposed as 

CG(r,l,n)/S n ~ V(k), 

ke.N(r,n) 

where V(k) is an irreducible representation of G(r,l,n) with dimV(k) = ( fc n k ) . 

(3) iyei 6e i/ze zona/ spherical function corresponding to V(k). For x = (xi, ■ ■ ■ ,x n ;a) € 
G(r, l,n), we have 

, , m A(fc)(?l'-- - >Zn) 
= 7 n \ ' 

Vfeo/-- ,fer-l' 
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where £j = £ x * . Moreover, the table of the zonal spherical functions is given by 

( m x(k)( £ ) \ 



( n ) 

\k ,-,k r -iJ / k,£eN(r,n) 



m » (£) 

We denote by lo^ £ = -, — — y the value of zonal spherical function indexed by (k,£). 

3. B(r,n) and (G(r,l,n), S n ) 

We consider a stochastic processes on B(r,n). The initial distribution is given by vo(b) = 

1, b = Jo, , 

Let 7r be the uniform distribution on B(r,n), i.e. it = — . Similarly we denote by 

0, 6^ J . 
tt the uniform distribution on G(r, 1, n). 

Let P = (p(a, &)) a ,6e_B(r,ra) be a G(r, 1, n)-invariant stochastic matrix on P(r, n) i.e., p(ga, 56) = 
p(a,b) for any a,b £ B(r,n) and x € G(r, l,n). Then we define a function ^ on B(r,n) by 
1^(6) = p(bo,b). Since the action of G(r, l,n) is transitive, there exists 5 G G(r, l,n) such that 
6 = g&o- Put = ^jp(bo,gbo). It is easy to check that ^ is a stochastic distribution and a bi-5 n 
invariant function on G(r, 1, ra). Therefore v can be expanded by the zonal spherical functions, 
say v = J2keN(m) a k UJ k- Then the orthogonality relation of the zonal spherical functions gives 
us the following proposition. 



Proposition 3.1. Let k = (k ,--- ,k r -i) G N(r,n). Put f^ = YlgeG(r,i,n) "(d) 00 ^ )■ ^ en 
i/ie coefficients 's are expressed by 



( n ) 

Vfc o,-" ,k r - l/ 

r n n\ 

Proof. Take an inner product of the both sides of v. □ 



The probability being in a state b = gbo (g G G(r, l,n)) after iV-steps iterate with a start 
point 60 is 

P {N) (b)= Yl p(b Q M)p{biM)---p{b N -i,b) 

bo,--- ,b k -ieB(r,n) 

= P{ b o,gibo)p{gib ,g 2 bo)---p{g N -ib ,gb) 

bo,--- ,b N - 1 eB(r,n) 

= ^2 v{xi)v{x^[ 1 X2) ■ ■ ■ ^(xj^^g) = v* N {g) (iVth. convolution power), 

XI,-" ,x n -i&G(r,l,n) 

where we denote by gi an element satisfying bi = gib®. By using a property of the zonal spherical 
functions Wfc * Wfe , = Mggg^,, we have 

( n ) 

~*N „N,*N Vfe °'"'" ' k r-i' fN, , 

" = a fc w fc = — — 4 w fc- 

keN(r,n) keN(r,n) 
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For stochastic distributions \i and // on a space X, the total variation distance is defined by 



Then the following estimate is known. 
Proposition 3.2. [1, Corollary 4.9.2, pp. 144] 

\\,,*N m2 1 1 ~*JV ~||2 



'lhv = ij>W-/<'M|. 



where k = (ko,--- ,k r -{) runs over N(r,n) except k = (n, 0, •■■ ,0) which corresponds to the 
trivial representation of G(r,l,n). 

We compute /^.'s for three distinct stochastic matrices on B(r,n). We define a function on 
B(r,n) x B(r,n) by 

dr(6,c) = #{i | 6j / Cf}, 
where 6 = (bi | 1 < i < n) and c = (q | 1 < i < n). We define a stochastic matrix by 



p(x,y) 



d r (x,y) = l 
0, d r (x,y)^l. 



Clearly the setting means the first way of shuffle explained in Section 1. 
Proposition 3.3. For k = (ko, • • • , fc r -i) € B(r, n), we have = ( £™ — 1 



j {i = m), 

Proof. We define g m j = (gi, ■ ■ ■ ,g n : 1) € G(r, l,n) by & = < for any 1 < m < n 



(i ^m) 



1 (9 = 9m,j), 



and 1 < J < r — 1. Then we have z/(p) = < ' Now we can compute f^ as follows. 

(g^g m j). 



r—l n 



Yl %H (9) = re(r 1 _ 1) Yl E U \(k) (9mj ) 
g€G{r,X,n) 3 = 1 ™ =1 

r-1 %r^r— 1 / n—1 \f~ i 3 -. 1 r r—l 

1 r^LHlJo,-^-!,-^.!] 1 ! 



(. " ) ~ r - 1 ^ n ^ ^ 

j = l Vfco,-" ,Kr-l' 1=0 j = l 

1 f fed (r-1) | y^-fe^l = 1 /rfcp \ 
r — l 1 n r-^ ra j r — 1 \ n / 



In the second equality we use Theorem 12.31 (3) and the following equation; 



□ 
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From the proposition above, we can easily show \f^\ < 1 and \f^\ = 1 O ko = < 
Theorem 3.4. For g = (x\, ■ ■ ■ ,x n : a), we have 



n (r > 3) 

or n (r = 2). 



fc =o v v y/ 

/ ,v I — 1 {xi = 1), 
w/iere <3?j = £ Xl + + • • • + ^ r >Xi = < and ej(xi, • • • ,x n ) is the j-th ele- 

mentary symmetric polynomial. 

Proof. For k = (ko, k\, ■ ■ ■ , k r -i), we have the following generating function ([3]) 



Ila +*«) = £( 

i=l k V 



ko, ■ ■ ■ , k r —i 

We remark that niL^i+^O = Ej e i(^ir-- , $n) ande, ($i, • ■ ■ , $ n ) = E^(A(fc))=j m A(fc)(£i> " ' ' 
where is an elementary symmetric polynomial and £(A) is the length of A. We compute 



1 U ,--- ,fc r _i) fN^j 

keN(r,n) 

r"7?J 2 (n(r-l)) 2 m A(fc)<&> ' ' ' ^») 

fc =0 v v y/ feiH \-k r -x=n-ko 

1 / rko - n \ N 

r n n! ^-^ V n r — 1 / 



fc =0 

□ 



From this theorem, we have immediately the following corollary. 

-*7V _ 1 

urn 1, then 



Corollary 3.5. If r > 2, then lmiAr^ 00 z> = If r = 2, let £ be a number of balls in the 



lim ~*2N = I M (* = (mod 2)) ^ ^-1 = J 2=^ST (/ = ! ( mo d 2)), 
N -">° 1o (£=1 (mod 2)) ' N ^°° lo (£ = (mod 2)) 

We try to estimate an upper bound of v* N . Before state a theorem, we see the following 
example. 

Example 3.6. We put r = 3 and n = 20. Then we have the following the graph of total 
variation distance. 
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2.0 



1.5 



1.0 



0.5 




Here The horizontal axis is the number of shuffles N . 
Theorem 3.7. Put N = (log r n + c) . 

(1) In the case o/r = 2, we have 

\W* N -A\ 2 TV-V±<\e- C - 

(2) In the case of r > 3 ; we have 

\\,,*N i|2 ^ 1 -c 

\\v -TT\\ TV <-e . 
Proof. We refer to Proposition 13.21 and compute 




Here we use |1 — a\ < e a for a < 1. □ 

Remark 3.8. In the Theorem above, "1/4" is the limit of \\u* — vr||yy as n — > oo which comes 
from a consequence of Corollary \3. 51 



4. Appendix 

Through this section, the urns Uq, U\, ■ ■ ■ , U r -i form a circle. Here we take up another shuffles 
and compute their Fourier coefficients. 
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4.1. cyclic shuffle 1. We consider the following type of shuffle; "The ball picked randomly is 
transferred into the next urn on the left." We define the following stochastic matrix on B(r,n); 



p(b, c) = < 



d r {b, c) = 1 and 6; - a G {0, 1, r - 1} (1 < Vi < n), 
0, otherwise. 



Clearly this gives us the probabilities of the shuffle introduced above. Now we have 
Proposition 4.1. For k = (fco, • • • , fc r -i) G B(r, n), we have = Y7i=o 7rC _ *- 

Proof. We define g mJ = (g u ■ ■ ■ , g n : 1) G G(r, 1, n) by ^ = jS im . 

I n 

Yl v(9)u k {g) = ~Y1 W A(fe)(5m,l) 
geG(r,l,n) m=l 

Er—1 I n—1 r—1 , 

Vfco,— ,k r -i> i=0 



□ 



4.2. cyclic shuffle 2. Here a shuffle, "The ball picked randomly is transferred into the next 
urn on the right or left randomly," is considered. The stochastic matrix on B(r, n) of this shuffle 
is given by 



p(b,c) 



±, d r {b, c) = 1 and b { - a G {0, ±1, ±(r - 1)} (1 < Vi < n), 
0, otherwise. 

Now we have 

Proposition 4.2. For k = (ko, • • • , &V-i) £ -B( r > w), we Ziawe = ^-r — 1 
Proof. We define g mJ = (g u ■ ■ ■ , g n : 1) G G(r, 1, n) by & = j'<S im . 

1 " \ 

%)Wfc(</) = ^Y1 ( W A(fc)(5m,l) + W A (fc)(ym,r-l) y 



9 eC7(r,l,n) m=l 

^r— 1 / n— 1 



l ELo (fc ,-,fc^i,-,fc r _i) (C * + = i^ , j ^ 

2 (. \ ) 2 ^ n l ^ + W " 

Vfco,— ,K r -l/ j=0 



□ 
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